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Abstract 



m 
o 

We have constructed star models consisting of four parts: (i) a homogeneous inner core with 



anisotropic pressure (ii) an infinitesimal thin shell separating the core and the envelope; (iii) an 



^ . envelope of inhomogeneous density and isotropic pressure; (iv) an infinitesimal thin shell matching 

the envelope boundary and the exterior Schwarzschild spacetime. We have analyzed all the energy 
conditions for the core, envelope and the two thin shells. We have found that, in order to have 
static solutions, at least one of the regions must be constituted by dark energy. The results show 
that there is no physical reason to have a superior limit for the mass of these objects but for the 
ratio of mass and radius. 



*Electronic address: chan@on.br 
^Electronic address: mfasnic@gmail.com 



^Electronic address: jaime@mast.br 



1 



I. INTRODUCTION 



Over the past decade, one of the most remarkable discoveries is that our universe is 
currently accelerating. This was first observed from high red shift supernova la [L-El anc ^ 
confirmed later by cross checks from the cosmic microwave background radiation 8j- 9] and 
large scale structure [h]]- 1o1 |. 

In Einstein's general relativity, in order to have such an acceleration, one needs to intro- 
duce a component to the matter distribution of the universe with a large negative pressure. 
This component is usually referred as dark energy. Astronomical observations indicate that 
our universe is flat and currently consists of approximately 2/3 dark energy and 1/3 dark 
matter. The nature of dark energy as well as dark matter is unknown, and many radi- 



cally different mode 



s have been proposed, such as, a tiny positive cosmological constant, 

re non-linear F(R) models 21J- 23|, and dark 



24j- 43]; see also the review articles 44]- 



45|, 



quintessence [16fl-|l8j, DGP branes [l9j-[20(, t 
energy in brane worlds, among many others 
and references therein. 

On the other hand, another very important issue in gravitational physics is black holes 
and their formation in our universe. Although it is generally believed that on scales much 
smaller than the horizon size the fluctuations of dark energy itself are unimportant 46], 



their effects on the evolution of matter overdensities may be significant [47|]- 48] . Then, a 
natural question is how dark energy affects the process of the gravitational collapse of a star. 
It is known that dark energy exerts a repulsive force on its surrounding, and this repulsive 
force may prevent the star from collapse. Indeed, there are speculations that a massive star 
doesn't simply collapse to form a black hole, instead, to the formation of stars that contain 
dark energy. In a recent work, Mazur and Mottola [49] have suggested a solution with a 
final configuration without neither singularities nor horizons, which they called "gravastar" 
(gravitational vacuum star). In this case, the gravastar is a system characterized by having 
a thin shell but not infinitesimal made of stiff matter, which separates an inner region 
with de Sitter spacetime from the Schwarzschild exterior spacetime. The elimination of the 
apparent horizon is done using suitable choice of the inner and outer radius of the thin shell, 
in such way that the inner radius be shorter than the horizon radius of de Sitter spacetime 
and the outer radius longer than the Schwarzschild horizon radius. In a later work, Visser 
and Wiltshire 50] have shown that the gravastar is dynamically stable. The possibility of 
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the existence of objects like the gravastar brings all the discussions about the fact that is 
unavoidable that gravitational collapse always forms a black hole. As a result, black holes 



may not exist at all 51 ]- 



521 ] . Another related issue is that how dark energy affects already- 



formed black hole is related to the fact that it was shown that the mass of a black hole 
decreases due to phantom energy accretion and tends to zero when the Big Rip approaches 



531 ]- 54|. Gravitational collapse and formation of black holes in the presence of dark energy 



were first considered by several works 



55j-[58|. 



Based on the discussions about the gravastar picture some authors have proposed alter- 
native models 5t|- 63]. Among them, we can find a Chaplygin dark star 64], a gravastar 
supported by non-linear electrodynamics 65j, a gravastar with continuous anisotropic pres- 



sure 



661 ] . Beside these ones, Lobo [67| has studied two models for a dark energy fluid. One of 



them describes a homogeneous energy density and another one which uses an ad-hoc mono- 
tonic decreasing energy density, both of them with anisotropic pressure. In order to match 
an exterior Schwarzschild spacetime he has introduced a thin shell between the interior and 
the exterior spacetimes. 

Since all the works cited previously are based in particular solutions, the investigation 
of more general solutions or others particular ones is important in order to establish the 
generality of these results. Our aim with this work is to construct another alternative model 
to black holes considering the possibility of gravitational trapping of dark energy by standard 
energy, i.e., not dark energy. In this context we could have the dark energy sustaining the 
collapse of the standard energy, while the standard matter would trap the dark energy. In 
our model the mass function is a natural consequence of the Einstein's field equations and 
the energy density as well as the pressure decreases with the radial coordinate (envelope), as 
expected for known stellar models. In order to eliminate the central singularity present in 
this model, we have considered a core with a homogeneous energy density, described by the 
Lobo's first solution 63]. The junction between the envelope and the Schwarzschild exterior 
spacetime has imposed a presence of a thin shell, as well as, the junction between the core 
and the envelope. 

The paper is organized as follows. In Section II we show the Einstein field equations. In 
Section III we present a particular solution that represents an isotropic and inhomogeneous 
dark energy fluid. In order to construct a more realistic model we consider it as the envelope 
of the star with a core constituted by a regular anisotropic and homogeneous fluid and a 
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Schwarzschild exterior. In Section IV we show the junction conditions between the core and 
envelope regions and between the envelope and the exterior spacetime. Finally, in Section 
V we present our final considerations. 



II. THE FIELD EQUATIONS 

We consider here a static spherically symmetric spacetime given by the following metric 

/dv 
g(r)dr dt 2 + — — -— + r 2 (d6 2 + sin 2 6d<f) 2 ), (1) 
J 1 — 2m{r)/r 

where g(r) and m(r) are arbitrary functions of the radial coordinate, r. 

The stress-energy tensor for an isotropic distribution of matter is given by 

T„ v = (p + p)v yVu + PSV> ( 2 ) 

where v M is the four-velocity, p is the energy density, p is the radial pressure measured in 
the radial direction. 

Thus, the Einstein field equation, G^ v = 8nT^ u , where is the Einstein tensor, provides 
the following relationships 

m! = 47ir 2 p, (3) 

m + 4nr 3 p , d . 

9 = —, 7, — r> ( 4 ) 

r[r — 2m) 

i (p + p)(m + 4irr 3 p) 

P = 7 7T-\ ' 5 

r[r — 2m) 

where the prime denotes a derivative with respect to the radial coordinate, r. Equation (J5]) 
corresponds to the Tolman-Oppenheimer-Volkoff (TOV) equation. 

The standard cosmology considers the Universe constituted of a perfect fluid described 
by an equation of state of the type p = up. Using this type of equation of state with a 
sufficiently negative pressure could explain the positive acceleration of the Universe in the 
context of the General Relativity theory. We consider this kind of equation of state in order 
to investigate the conditions of existence of dark energy star. Assuming an equation of state 
of the type p = up and substituting it in the equation (jSJ) we have 

v+( P+ ^) (m . +4 ; r3 T ) =°- (6) 

r[r — 2m) 
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Solving equation (jSJ) in terms of m we get 

wr 2 (p' + A-nrujp 2 + Anrp 2 ) , 

m = . (7) 

2u;p / r — up — p 

Differentiating equation © and substituting into equation (j3J) we obtain 

— f-247rr 2 u; 2 p 2 p / - 12-Kr 2 p 2 up' + 287rru; 2 p 3 + 207rrcjp 3 + Arcrp 3 

(2ujp'r — uop — p) 1 v 

- 3u 2 p' 2 r + 2p'cj 2 p + 2p up - cup ,2 r + f 27rruJ 3 p 3 + rp"u 2 p + %ixr 3 u 3 p 2 p" + %ixr 3 u 2 p 2 p" 

- 167cr 3 u 3 pp' 2 - f 2-Kr 2 uo 3 p 2 p' + rp"wp) = 0. (8) 

Below we work with a particular solution of this equation. 

III. SOLUTION OF THE PHYSICAL EQUATIONS 

Solving equation (jSj) in terms of p we can obtain a particular solution written as 



to 

2n(co 2 + 6u + l)r 2 



P= «-t.x , ... ■ ( 9 ) 



where we can note that the values of a? can not be — 3 + 2^2 or -3 - 2-^2, in order to 
avoid the singularities for these values. Since p > thus —3 — 2\/2 < < —3 + 2y / 2 or 
a; > 0. Besides, our envelope with u — 1 reproduces the same rigid fluid used by Mazur and 
Mottola [49], with the pressure proportional to r~ 2 . 
Thus, using equations (j3J) and (jlj) we obtain that 

^ l^TTv < 10 > 

The metric can be rewritten as 

ds l = _ r ^+D dt 2 + u}2 + 6uJ + 1 d r 2 + r 2 (d6 2 + sin 2 6d<p 2 ) , (12) 

(c<j + l) 2 v 7 

where the subscript (+) denotes the envelope spacetime as we consider it (see below). Since 

the metric function g rr is negative for —3 — 2y/2 < cu < —3 + 2y / 2, this interval is not 

allowed. Thus, this solution implies uo > 0. 

The Kretschmann scalar for ffl2l is given by 

K-R Jr#r - ^W + 2^ + 7) 

K ~ Ra ^ R ~ r^ 2 + 6u + l) ' (13) 



where Rappv is the Riemann tensor. We can note that lim^oo K = 0, i.e., this solution is 
asymptotically flat. But, we can see also that Um r ^,oK ~~ > +00, the solution is singular at 
the origin. Then, we consider the metric (I12p as an envelope solution. Below we present a 
core solution. 



A. Core Solution 

In order to avoid the singularity at r = we cut the spacetime (fT2|) around its origin and 
fit another one, an anisotropic fluid with the density of energy u constant. We have chosen 
an anisotropic solution for the core since it was shown that gravastars solutions must exhibit 



anisotropic pressures to be finite-sized objects 



661 ] . Using the results from Ref . [67| we have 



and 



Pt = Pr 



H — l^o = constant, 
p r = kfj, Q , 

~ (l + A;)(l + 3A;)M(r) 



(14) 
(15) 



(16) 



2k(r - 2M(r)) 

In this case, we have as isotropic pressure limits, k — — 1/3 and k = —1, otherwise we have 
anisotropic pressures. Here we point out that the anisotropy in the pressure changes the 
ranges where the energy conditions are satisfied, which depend on the values of M(r)/r [68]. 
The core metric can be written as 



ds< 



2M(r) 



-1 -(l+3fc)/2 



dv 2 



2M(r) 



1 -1 



dr 2 + r 2 {d0 2 + sin 2 Odcf) 2 ), (17) 



where M(r) = 47r/i r 3 /3 and the subscript (-) means the core spacetime. See figured! 

Note from the equation ( fl6|) that the anisotropy of the core \pt — p r \ decreases with the 
r coordinate. Then, the supposition of an isotropic envelope is physically reasonable. 

B. Exterior Solution 



In order to limit the matter of the star we match the exterior spacetime with a 
Schwarzschild solution. The Schwarzschild exterior metric can be written as 



dsi 



dr 2 + f 2 {d6 2 + sin 2 



(18) 



EXTERIOR 




EXTERIOR 



FIG. 1: The regions of the star: the core, the envelope and the exterior. The radius of the core is 
rs and the total mass inside i*e is M(rs). The density of the core is homogeneous and the pressure 
can be isotropic or anisotropic. The outer radius of the envelope is Fg. 
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where M is the total mass of the star, including the mass of the core, and the subscript (e) 
denotes the exterior spacetime. 

IV. JUNCTION CONDITIONS 



The metric of the hypersurface £ at the frontier of the core and the envelope is given by 



dsl = -dr 2 + R(T) 2 (d6 2 + sin 2 6d<j) 2 ), 



(19) 



where r is the time coordinate define only on S. 

The metric of the hypersurface £ at the frontier of the envelope and the exterior is given 

by 

ds\ = - df2 + R{r) 2 {d6 2 + sin 2 6d<f) 2 ), (20) 

where f is the time coordinate define only on £. 

From the conditions (ds 2 _ ) E = (<is+) E and (ds 2 + )f, = (ds 2 ) E and using equations f|T2l) . 
(|T7|) and (TLB]) , we get the following relations 



r E = r E = R(t), 
f E = r E = R(f), 



dv s 
dr 



1 2M(r s 
r E 

\ _ -4u/(w+l) 



(l+3fc)/2 



df 



1 - 



2M 



(§)- 

The core extrinsic curvature is given by 



r S / 

4w/(tj+l) 



(21) 
(22) 
(23) 

(24) 

(25) 

(26) 



2r s 



-2M(r E ) + r E 
r E 



-i -3fc/2 



dr 2M(r E )-r E 



[M(r E ) - M'(r E )r E ] (3fc + 1), 



-2M(r E ) + r E 



-1/2 



[2M(r E ) - r E ] 



k: 



-2M(r E ) + r E 



r E 



[2M(r E ) - r E ] sin 2 0, 



(27) 
(28) 

(29) 
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and the envelope extrinsic curvature is given by 



K+ = -2r- 



-4/(w+l) 3 



dt 



Vou 2 + 6u + l\ — 



dr ) uj 2 + 6u + 1 
1 



uj 2 + 6cj + 1 
1 



a; 2 + Quo + 1 



rs sin 9. 



The core extrinsic curvature can be rewritten as 

(1 + 3fc) " 



k: 



2r| 



l 2M(r s ) 



= -<l + 3*)^ 



1 



K 0B = r S 



1 - 



[M(r E ) - M'(r s )r s ] 
2M(r s ) ' 

2M(r s )~ 



^ = ^ sin 2 9, 
the envelope extrinsic curvature can be rewritten as 



= r^{uo + 1) (uj 2 + Quo + , 
ff+ = tf+ sin 2 61, 

and the Schwarzschild exterior extrinsic curvature can be written as 



'T -2 



K e ee = f fi 1 

K*=KLsm z 9. 



2M 



2M 



,2 , 



A. Junction conditions between Core-Envelope: without a thin shell 



Using the junction conditions K TT = K+ T and K ge = K^ e we get two equations in terms 
of M(r s ) and r s 



-(1 + 3*)^ 



2M(r s ) 



H (u 2 + + 1 



0. 



and 



1 - 



2M(r s ) 



r s (o; + l) (cj 2 + 6cj + 1 



0. 



(42) 



(43) 



We can solve these two equations, obtaining the total mass of the core and the radius of 
the core as a function of uj and k giving 



us = 3k, 



(44) 



and 



M(r E ) 



2u 



< 



(45) 



r E (u 2 + Qu + 1) ~ 4' 
Note that although there is a negative interval for u where the mass is positive, this same 
interval is forbidden by the metric's signature, equation (Tl2l . Consequently, k is also non 
negative. We note that the solution can represent a star model with an upper limit for the 
ratio of mass to radius of the core. 



B. Junction conditions between Core-Envelope: with a thin shell 



So, we match a thin shell of energy density a and pressure P at the frontier between the 
core and the envelope. Thus, we can write 

1 a 



where 



Since 



a = Kn 

4?r 8 



r s v ' rf; 



\ 2M(r s ) 
rs 



(46) 
(47) 

(48) 
(49) 
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then 



P 



87rr E 



+ 1) (u 2 + 6u + 1 



1 



2M(r s ) 



rs 



(7 



4vrr E 



1) (u 2 + 6uj+ l)' 



1 - 



2M(r s ; 



(3w + 1) (w 2 + 6cj + 1 



1 + (1 + 3A;) 



M(r s ) 



2M(r s ) 
rs 



(50) 



(51) 



(52) 



where < — 3 — 2v^2 or u; > — 3 + 2y/2 (in order to have a and P real) and r s > 2M(rs). 
Since a > and p > then u; > 0, which represents a not dark energy envelope, but no 
restriction is imposed to the core, i.e., we can have any values for k. 



1. Energy conditions for the thin shell: Core-Envelope 



The energy conditions 69] for a thin shell can be written as 

1. Weak: a > and ff + P > 

2. Dominant: a + P > and a - P > 

3. Strong: a + P > and a + 2P > 

The characterization of dark energy fluid is the violation of one of the strong energy 
conditions, more specifically, that one related the Raychaudhuri equation 68( . If the second 
of the weak energy conditions is violated, we have a phantom dark energy fluid. 

For the inner thin shell, that is, between the core and the envelope, we have 



a + P 



o-P 



87rr s 
1 

8vrr s 



10J 



6u 



- 1) 

(huj + 3) (uj 2 + Geo + 1 



a > 



+ 



(5 + 3k) 



M(r s ) 



rs 



2M(r s ) 



cr + 2P 



1 



47rr s 



2u u 2 + Qu + 1) -3 



k) 



(1-k) 



M(r s ; 
rs 



M(r s ) 
rs 



rs 
2M(r s ) 



(53) 
(54) 



rs 



1 - 



2M(r s ) 
rs 



(55) 
(56) 



In the Table I, we summarize the results of the energy conditions of the inner thin shell for 
different limits and k > —1/3. We analyzed only the case where k > -1/3 because we are 
interested to check if it is possible to have all the structures constituted by standard energy. 
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TABLE I: This table summarizes the results of the energy conditions of the inner thin shell for 
different limits and k > —1/3. 



Case 



87rr s ((T + P) 



8vrr s (cj - P) 



4vrr E (o- + 2P) 



uj -> 
M(r s )/r s < 1 



.(4 + 3A;)M|^ < o 



3A;M(££) > o 



. 3 (i + A; )M^l<0 



w -> 



M(r s )/r s 1/2 



4 l-(5+3fc)/2 < q 

/ 1 2 M(r s ) - 



o 3[l-(l-fc)/2] > Q 



1-2 



M(r a ) 



-3(l + fc)/2 



A/(r S ) 



< 



uj -> 1 



M(r s )/r E < 1 



l-(4 + 3£;)^g^ >0 

M(r S ) <■ 1 
r E — 4+3/c 



V2 



+ 3 



1 + fc 



M(r s ) 



>0 -^-3(l + fc)^i>0 



M(r s ) ^ J_ /_4_ _ o 
rs ~ 3fc 



M(r s ) 



< 



3>/2(l+fc) 



W -> 1 



M(r s )/r s -► 1/2 



l-(5+3fc)/2 

1 2 A//(r s )" 



< 



_J_ 3[l-(l-fc)/2] > Q 

a/2 / 1 2 M(r s ) - 



J_ _ 3(l+fc)/2] 
V2 



1-2 



M(r S ) 



< 



C. Junction conditions between Envelope-Exterior: without thin shell 

Using the junction conditions K^- = K~ and Kg~ 6 = Kg B we get two equations in terms 
of M and r s 



M 



and 



1 - 



2M 



2M 



2 2cj / 9 
+ — [uj 2 + 6cu + 1 



0. 



2 - Vyiuj + 1) f^ 2 + QUO + 1^ 7 



0. 



(57) 



(58) 



Considering these two equations we have the equation 



(59) 



which gives us 



u = 0, 
= -2, 



(60) 



where only the first one is solution of the original equations (1571) and (1581) . This solution 
implies that 

M = M(r s ) = 0, (61) 
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which means that all the spacetime is Minkowski. Then we can conclude that for this kind 
of fluid distribution, it is not possible to have structures surrounded by a Schwarzschild 
spacetime without the introduction of a thin shell. 

However, in a Lobo's work [67| it was necessary the introduction of a thin shell on the 
junction hypersurface. The authors have suggested that it is possible to match the interior 
and exterior spacetime without a thin shell. However, eliminating the thin shell, vanishing 
its energy density and pressure (equations 25 and 26 in the original work), we get the same 
junction conditions considered in this present work, i.e. giving a Minkowski spacetime. In 
order to do that we assume m = M and u = or m' = with also a = 0. Thus, the 
unique interior solutions that admits a matching with Schwarzschild spacetime, without the 
introduction of a thin shell, is the Minkowski (wf = 0) and the dust (ou = 0) solution. 



D. Junction Conditions between Envelope-Exterior: with a thin shell 



Since as seen in the above subsections, we can not match an exterior Schwarzschild 
spacetime directly with the envelope, we build a thin shell of energy density a and pressure 
P in order to do this match. Thus, 

1 a 



where 



a 



"to** 



1 



K e - K- -. 

13 v 



(62) 
(63) 

(64) 



Since 



then 



- M ( 
kI = -7T 1 - 



a 



2M\ 2 2u , ., 
— -— uj z + Qlu + 1 



2M\ 2 1 . . / „ 

— -— (uj + I) (u 2 + Qu + I 



P 



2M\ 



1 r s) V 1 ^ J 



10J 



+ 1) (u 2 + 6u + 1 



- (3a; + 1) (u 2 + 6^ + 1 



(65) 
(66) 

(67) 
(68) 
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where uj < —3 — 2y2 or uj > —3 + 2 a/2 (in order to have a and P real) and f g > 2M. Since 
p > then the unique physical system must have uj > 0. 
Since a > then we must have 

2uj 



for cj > 0, and where 



M 1 
+ 6uj + 1 ~ ?s" ~ 2' 

2cj 



1 

- < -. 

uj 2 + Quj + 1 ~ 4 



(69) 



(70) 



Energy conditions for the thin shell: Envelope-Exterior 



The energy conditions 69j for the thin shell between the envelope and the exterior can 
be written as 

a > 

1 



a + P 



a 



- P 



1 



57rri 



-(5w + 3) (uj 2 + Quj+1) 5 + 
(u - 1) (cj 2 + Quj + 1 





3-5*" 




2M~ 




+ 






1 






r 2_ 




r 2 . 





(71) 
(72) 



+ 



a + 2P 



1 



AlXT% 



_ 1 

-2(2cj + 1) (cj 2 + 6c; + l) 2 + 



M" 




2M" 




1 - 3— 




1 


") 


r 2_ 




r 2 . 




M 




2M 




2-3— 




1 




r 2. 




r 2 . 





(73) 



(74) 



In Table II we summarize the results of the energy conditions of the outer thin shell for 
different limits. 

Comparing the Tables I and II we can have the following conclusions: 

1. Limits uj -> 0, < 1, ¥■ < 1: both (inner and outer) thin shells are phantom; 



2. Limits a; — > 0, 



M(r s ) 



1/2, ^ — > 1/2: the inner thin shell is made of dark energy 



and the outer thin shell violates the dominant energy condition; 
3. Limits uj -> 1, < 1, # < 1: the inner thin shell is constituted by dark and 



not dark energy, depending on and |f , and the outer thin shell is made of dark 



energy; 



4. Limits a; — ► 1, 



M(r s ) 



1/2, — > 1/2: the inner thin shell is made of dark energy 



and the outer thin shell is constituted by not dark energy, violating the dominant 
energy condition. 
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TABLE II: This table summarizes the results of the energy conditions of the outer thin shell for 
different limits. 



Case 


87rr s (a + P) 


8?rr s (CT - P) 


4vrr s (a + 2P) 


uj -> 


-2M- < 


-2*L < 


-#<o 




r E 


r E 


r E 


M/fg < 1 









uj^O -3 + ^2: — >0 -1 JJl — <0 -2 iZ2 — >0 

1-2 M- J1-2M- J1-2M- 

r E V r S V r S 



M/fg -» 1/2 



M/rg < 1 



M/rg -> 1/2 



r s - 2V 7 ! rg - 2 v 7 ! r E ~ 



V2 /i_ 2 jy. Jl-2M- ^2 J1-2M- 



E V E V E 



We can conclude with these limits that we always have one or both thin shells constituted 
by dark energy. 

V. CONCLUSIONS 

We have constructed a star model consisting of four parts: (i) a homogeneous inner core 
with anisotropic pressure (ii) an infinitesimal thin shell separating the core and the envelope; 
(hi) an envelope of inhomogeneous density and isotropic pressure; (iv) an infinitesimal thin 
shell matching the envelope boundary and the exterior Schwarzschild spacetime. We have 
analyzed all the energy conditions for the core, envelope and the two thin shells. 

In our model the mass function is a natural consequence of the Einstein's field equations 
and the energy density as well as the pressure decreases with the radial coordinate (envelope), 
as expected for known stellar models. In order to eliminate the central singularity present 
in this model, we have considered a core with a homogeneous energy density, described by 
the Lobo's first solution [otI ]. 

We have proposed in this work an alternative model and a generalization to gravastars. 
Note that for k = — 1 we get a vacuum core with a cosmological constant A = 87Tfio, i.e., a 
de Sitter solution. Thus, we have constructed models with the same structures of the Mazur 
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and Mottola's one [491 ]. with five regions (an interior de Sitter solution + an infinitesimal 
shell + a non infinitesimal shell + infinitesimal shell + an exterior Schwarzschild solution). 
However, in our models each one of these regions is made of more general kind of fluids. 

In the Table III we summarize the results of the match of the core's spacetime with 
envelope's spacetime and the results of the match of the envelope's spacetime with the 
exterior spacetime. 

Combining the results of the Tables I, II and III we can see that in the analyzed cases we 
always have the presence of the dark energy at least in one of the thin shell or in the core. 
Note also that from the results that there is no physical reason to have a superior limit for 
the mass of these objects but for the ratio of mass and radius, in order to find out which 
one is made of dark energy. 
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